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Abstract. The percolation phase transitions of two-dimensional lattice networks under a generalized 
Achlioptas process (GAP) are investigated. During the GAP, two edges are chosen randomly from the 
lattice and the edge with minimum product of the two connecting cluster sizes is taken as the next occu- 
pied bond with a probability p. At p = 0.5, the GAP becomes the random growth model and leads to the 
minority product rule at p = 1. Using the finite-size scaling analysis, we find that the percolation phase 
transitions of these systems with 0.5 < p < 1 are always continuous and their critical exponents depend 
on p. Therefore, the universality class of the critical phenomena in two-dimensional lattice networks under 
the GAP is related to the probability parameter p in addition. 
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1 Introduction 

The percolation phase transition concerns the formation of 
a macroscopic component in systems on both lattices and 
networks [T] . It provides a model for the onset of a macro- 
scopic component in random media T and social networks 
[5]. It was widely believed that the percolation transition 
is a typical continuous phase transition for various net- 
works [3]. However, Achlioptas, D'Souza, and Spencer [i] 
found recently that the percolation phase transition in 
random network becomes discontinuous (first-order) un- 
der the Achlioptas process (AP), where the edge with min- 
imum product of cluster masses is connected from the two 
randomly chosen unoccupied edges. The Achlioptas pro- 
cess suppresses the appearance of larger cluster and dis- 
courages the formation of a giant component, which has 
the size comparable with the number of vertices N . The 
percolation phase transition is delayed by this Achlioptas 
process and becomes sharper. It was argued by themf4] 
that the phase transition is discontinuous and named as an 
explosive percolation. Later, the Achlioptas process was 
introduced to the two-dimensional regular lattice networks 
[S] and scale- free networks (Sl[7]. It was claimed that the 
explosive percolation was found both in lattice [5] and in 
scale- free networks [61I3- 

In the article of Achlioptas et al.[4 , the step interval 
A between the size of the largest component Si = N-^^^ 
and Si = 0.5N is used as the criterion of continuous or 
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discontinuous phase transition. Later, Ziff applied this cri- 
terion to the two-dimensional regular lattice networks [SJ 
|8]. It was found that the size-dependence of A in lattice 
network is quite different from that in random network. 
It is not well established that the first-order phase transi- 
tion can be distinguished from the continuous phase tran- 
sition by the size-dependence of A. It was argued by da 
Costa et al.[5] that the explosive percolation transition, 
under their modified Achlioptas process, is actually con- 
tinuous. Recently, Riordan et al.piJj show mathematically 
that all Achlioptas process have continuous phase transi- 
tions. The finite-size behavior of the order-parameter dis- 
tribution function has been used as the evidence of both 
discontinuous [11 and continuous |1 2) phase transition. It 
is also argued with the finite-size scaling that the explosive 
percolation is continuous P^fH] . This controversy about 
the character of explosive percolation is going on and calls 
for further investigations. 

In this paper, we investigate the percolation phase 
transition in two-dimensional lattice network under a gen- 
eralized Achlioptas process (GAP), which will be intro- 
duced in the next section. The generalized Achlioptas pro- 
cess is characterized by a probability parameter p. The 
GAP becomes the random growth model at p — 1/2 and 
the minority product rule at p — 1. Using the finite- 
size scaling analysis, our Monte Carlo simulation results 
demonstrate clearly that the percolation phase transition 
in two-dimensional lattice networks under the GAP is con- 
tinuous. It will be shown that the critical exponents and 
therefore the universality class of the continuous percola- 
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tion phase transition depend on the probabihty parameter 

V- 

Our paper is organized as follows. In the next sec- 
tion, we introduce a generalized Achlioptas process in two- 
dimensional lattice network. In Section 3, we investigate 
the critical points of two-dimensional lattice network un- 
der the GAP and their critical exponents with the use of 
finite-size scaling. In Section 4, the finite-size scaling func- 
tion of the ratio S2/S1 is obtained at different probability 
parameter p, where 6*2 and 5*1 are the size of the second 
largest and the largest cluster in the network. The univer- 
sality class of the critical points in our model is discussed 
in Section 5. Finally we make some conclusions in Section 
6. 
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2 Two-dimensional lattice network under the 
GAP 

We consider a two-dimensional square lattice with size 
L X L and periodic boundary conditions in both direc- 
tions. There are N — V^ vertices in this lattice. We in- 
troduce a generalized Achlioptas process for adding edges 
into this lattice. In the generalized Achlioptas process, two 
edges are picked up randomly at each step. Each edge is 
connected with two clusters. The edge with the minimum 
product of the cluster sizes is chosen and added into the 
lattice with a probability p, where < p < 1. Correspond- 
ingly, another edge is chosen with a probability 1 — p. At 
p = 0.5, the GAP is equivalent to the classic Erdos-Renyi 
(ER) rule, where edges are picked up randomly. The two- 
dimensional square lattice with the ER rule is actually the 
two-dimensional bond percolation (BP) model. The GAP 
at p = 1 is the product rule (PR) of Ref. [3] and our model 
becomes the PR model on the two-dimensional regular 
lattice. 

In our Monte Carlo simulations, there are N isolated 
vertices in a two-dimensional lattice at the beginning and 
then edges are added into the lattice through the GAP. 
With the edges added, we obtain a network in the lattice. 
The lattice network can be characterized by a reduced 
edge number r = Nr/N, where Nr is the number of the 
edges added. 

We have used the algorithm of Newmann and ZifF [151 
[TB] in our Monte Carlo simulations. For the investigations 
related only to the largest cluster of lattice networks in 
Figs. 3, 5 and 7, the linear sizes L = 32, 64, 128, 256, 512, 
and 1024 are taken. When the second largest cluster in 
the lattice is taken into account in addition, only three 
linear sizes L = 64, 128, and 256 are taken in Figs. 2, 
4 and 6. To get enough samples for the average of each 
simulation, different steps are taken for the simulation of 
different system size. In our Monte carlo simulations, we 
run 10, 000, 000 steps for L = 32 until to 6, 400, 000 steps 
for L = 1024. 

For the cluster ranked R and with size Sji{r, L;p), we 
defined its reduced size as 



Fig. 1. Reduced size of the largest cluster, defined in Eq. ^, as 
a function of the reduced edge number r at the lattice size L = 
1024 and probability parameters p = 0.5, 0.6, 0.7, 0.8, 0.9, 
and 1.0. 

In Figlll we shown the reduced size si(r, L;p) of the largest 
cluster. 

In the Monte Carlo simulations of these data, we take 
the lattice size L = 1024 and the probability parameter 
p = 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. At small r, the re- 
duced size of the largest cluster is nearly zero. When r is 
large enough, the reduced size si becomes finite. This in- 
dicates the formation of a macroscopic component. There- 
fore, there is a percolation phase transition in the lattice 
network. The transition value of the reduced edge number 
re depends on the probability parameter p. It is shown in 
FigHlthat Tc increases with the probability parameter p. 
This is plausible since a larger p means stronger suppres- 
sion of larger cluster and therefore the later appearance 
of a macroscopic component. In the following, we will try 
to verify that these percolation phase transitions are con- 
tinuous or not. 



3 Critical points of two-dimensional lattice 
networks under GAP 

If the percolation phase transitions above were continuous, 
the reduced size of the cluster ranked R should follow the 
finite-size scaling form pTlfTB] 



SR{r,L-p)=L-PI''~SR{tL^I'^-p), 



(2) 



where t = {r — rc)/Tc characterizes the deviation from the 
critical point fc and v is the critical exponent of the cor- 
relation length ^ — ^ol^l "^^ This finite-size scaling form 
is supposed to be valid in the asymptotic critical region 
where L ;» lattice spacing and |i| <C 1. Outside the asymp- 
totic region, additional correction terms should be taken 
into account. 

For the largest cluster of lattice network, we have the 
finite-size scaling form 



SR{r,L;p) = SR{r,L;p)/N. 



(1) 



si{r,L;p) = L-P'-h{tL^'''\p). 



(3) 
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In the bulk limit L ^ oo, the reduced size of the largest 
cluster becomes 

si(r, oo;p)=0 (4) 

for r < Tc and 



si(r, oo;p) = Qp t^ 



(5) 



for r > re- The emergent macroscopic component is char- 
acterized by the critical exponent /3. The smaller the crit- 
ical exponent j3 is, the larger is the macroscopic compo- 
nent. 

Near the critical point, the reduced size of the second 
largest cluster can be written also in a finite-size scaling 
form 

S2{r,L-p)^L-^'''UtL^''';p)- (6) 

Using Eqs. ([TJ ([3]), and ([6]), we can obtain the finite-size 
scaling form of the ratio 



At the critical point t = 0, the ratio 



^2/^1 



C/(0;p), 



(7) 



(8) 



which is independent of the system size L. Therefore, the 
curves of S2/S1 at different system size L have a cross- 
point sX r ~ Tc. The critical point corresponds the fixed 
point of S2/S1, which can be used to determine the critical 
point of our system. 

The logarithm of Eq. ([3]) can be expressed as 

lnsi{r,L;p) = -{/3/iy)lnL + \nSi{tL^^'';p). (9) 

At the critical point r = re, we have 



lnsi(rc,L;p) = -(/3/j/) InL + lnsi(0;p). 



(10) 



which is a straight line with respect to InL. We can use 
this property to determine the critical point r^ of the lat- 
tice networks also. From the slope of this straight line, the 
critical exponent ratio l3/v can be determined. 

In the following, we use both the fixed point of S2/S1 
and the linear dependence of In si on In L as the criterion 
to determine the critical point of the two-dimensional lat- 
tice networks under the generalized Achlioptas process. If 
we have reached the asymptotic critical region, both si 
and S2 satisfy the finite-size scaling form in Eqs. ([3]) and 
([6|). The critical reduced edge numbers Tc obtained from 
82/ Si and Insi should be equal. On the other hand, the 
consistence of Vc obtained from two different methods can 
be used as the indicator of the accuracy of our simulation 
results. 

For the generalized Achlioptas process with probabil- 
ity parameter p — 0.5, edges are added randomly into a 
two-dimensional lattice and our model becomes the so- 
called bond percolation model. It is well known that the 
bond percolation model has a continuous phase transi- 
tion. We can determine the critical point of this model 
by the two methods described above. In Figl2j the ratio 
S2/S1 is shown as a function of the reduced edge num- 
ber r at different system size L. A fixed point between 
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Fig. 2. Size ratio S2/S1 of the second largest to the largest 
cluster at p = 0.5 and different L. There is a fixed point at 
re = 0.5000 ± 0.0002. 
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Fig. 3. Log-log plot of the reduced size si of the largest cluster 
at p = 0.5 and different r. We have chosen linear sizes L = 
32, 64, 128, 256, 512, and 1024. 



r ~ 0.4998 and r — 0.5002 is found and is in full agree- 
ment with Tc — 1/2 of the bond percolation model. In 
Tabledl we denote the critical point obtained from S2/S1 
asri^^ =0.5000 ±0.0002. 

As we have discussed above, a critical point can be de- 
termined alternatively by the linear relationship between 
Insi and InL. In Fig|3l lnsi(r, L;p) at different reduced 
edge numbers are shown. At r = 0.4996 and r — 0.5004, 
the curves of lnsi(r, L;p) are curved and their curvatures 
have different sign. At r = 0.5, the curve of lnsi(r, L;p) 
becomes a straight line. Therefore, we obtain the criti- 



cal point r, 



(2) 



with 



.(1) 



0.5000 ± 0.0004, which is in agreement 
From the slope of lnsi(r, i;p) at Vc, we get 
the critical exponent ratio /3/i/ = 0.108. Our Monte Carlo 
simulation results agree very well with the exact results 
Tc = 1/2 and P/v — 5/48 of two-dimensional bond perco- 
lation model, for a review, see Ref. [H]. 
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Fig. 4. Ratio S2/S1 at p = 0.8 and difTerent L. There is a Fig. 5. Log-log plot of the reduced size si of the largest cluster 
fixed point at r = 0.5208 ± O.OOOf . at p = 0.8 and different r. We have chosen linear sizes L = 

32, 64, 128, 256, 512, and 1024. 



For probability parameter p = 0.8, the edge that mini- 
mizes the product of two connecting cluster sizes is added 
into the lattice with a probability 0.8 from two randomly 
chosen edges. The connection of smaller clusters is fa- 
vored. The critical point of this system is investigated 
by the fixed point of S2/S1 and the linear dependence 
of Insi on InL. In FigHJ the ratio S2/S1 is plotted as a 
function of the reduced edge number r at different sys- 
tem sizes L. A fixed point of S2/S1 is found. It is be- 
tween r = 0.5207 and r = 0.5209. So there is a continuous 
phase transition in this system and the critical point is at 



0.5208±0.0001. Alternatively, this critical point can 



be determined from lnsi(r, L;p). In FigO lnsi(r, L;p) at 
r = 0.5205, 0.5207, and 0.5209 are shown. The curvature 
of \nsi{r,L;p) is negative at r = 0.5205 and positive at 
r = 0.5209. The function In si (r, L; p) at r = 0.5207 can be 
described quite well by a straight line with zero curvature. 
So the critical point rj; = 0.5207 ± 0.0002, in agreement 
with Tc . The slope of the straight line at re = 0.5207 
gives the critical exponent ratio f3/i' — 0.081. 

At p = 1.0, our model becomes the two-dimensional 
lattice network under the Achlioptas process. The ratio 
S2/S1 of this model is shown in FiglB] for different sys- 
tem size L. There is a cross-point between r = 0.5265 and 
r = 0.5267, which corresponds to the critical point of this 
system. Therefore its critical point is at Tc = 0.5266 ± 
0.0001. The curves of lnsi(r, L;p) at r = 0.52651,0.52655, 
and 0.52659 are shown in Figl?! Their curvatures change 
with r from negative to positive. The function becomes 
linear with respect InL around r — 0.52655. So we get the 
critical reduced edge number rj; — 0.52655 ± 0.00004, 
which agrees with re given above. From the slope of 
In si (re, L;p) with respect to InL, the critical exponent 
ratio P/v — 0.064 is obtained. Our results of the critical 
point and the critical exponent ratio are in full agreement 
with the results of Refs.[5ji5jX- 

In Fig|51 we summarize the critical reduced edge num- 
bers re of two-dimensional lattice networks under the GAP 




0.5256 0.5262 0.5268 0.5274 0.5280 
r 

Fig. 6. Ratio S2/S1 at p = 1.0 and different L . There is a 
fixed point at r = 0.5266 ± 0.0001. 
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Fig. 8. Dependence of the critical reduced edge number re on 
the probability parameter p. The error bars of Monte Carlo 
simulation data are smaller than the symbol. 
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Fig. 10. Finite-size scaling function U{tL^''';p) of the ratio 
S2/S1 at p = 0.5. The critical exponent of two-dimensional 
bond percolation v = 4/3 is taken for the scaling variable. 
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Fig. 9. Dependence of the critical exponent ratio P/v on the 
probability parameter p. 



with different p. It is found that re increases with p. At a 
larger p, the continuous percolation phase transition ap- 
pears at a larger critical reduced edge number re and the 
formation of a giant component is delayed. 

In Figini the dependence of the critical exponent ra- 
tio /3/1' on p is shown. With the increase of p, the ratio 
P/v decreases. We will show in the next section that l/iy 
increases with p. So it can be concluded that the criti- 
cal exponent (3 decreases with p. Smaller /3 indicates the 
stronger emergence of a giant component in the networks 
after the percolation transition. Therefore, the increase of 
p results in the delayed appearance of a continuous perco- 
lation phase transition and the formation of a larger giant 
component at the same time. 



is near the critical point re- For a given probability param- 
eter p, the different curves of S2/S1 at different L collapse 
into a finite-size scaling function after using the scaling 
variable tL^''^, where 1/ is the critical exponent of cor- 
relation length. In the following, we will investigate the 
finite-size scaling function of S2/S1 for different p. 

At p — 0.5, we use the critical exponent of two-dimensional 
bond percolation model ly = 4/3 [IH] for the finite-size 
scaling function of S2/S1. After defining the scaling vari- 
able ti^/" with this value of v, our Monte Carlo simulation 
results at L = 64, 126, 256 collapse and we get the finite- 
size scaling function of S2/S1, which is shown in FigfTOl 

At p = 1.0, the critical exponent v is unknown. Ac- 
cording to the finite-size scaling form in Eq. ([7]) , the curves 
of 5*2/ 5*1 at different system sizes can collapse only when 
the correct critical exponent v is used for the scaling vari- 
able. This property can be used also for determining the 
critical exponent v. At 1/z/ = 0.93, the curves of S2/S1 
at L = 64, 128, 256 collapse into its finite-size scaling 
function, which is shown in Fig llll 

In FiglT2l we demonstrate the variation of the finite- 
size scaling function of S2/ Si with the probability parame- 
ter p. In the region before the percolation phase transition, 
the finite-size scaling function of S2/S1 increases with p. 
The second largest cluster in this region is more important 
at larger p. In the region after the percolation phase tran- 
sition, the finite-size scaling function of S2/S1 decreases 
with p. The largest cluster in this region is more domi- 
nant at larger p. To get the finite-size scaling functions 
at different p, the corresponding exponents of correlation 
length are determined and presented in Fig]13l and Tabel 

EI 



4 Finite-size scaling functions of S2/S1 

In the last section, we have mentioned that the size ratio 
S2/S1 follows the finite-size scaling form in Eq.© when r 



5 Universality classes 

The concept of universality plays a fundamental role in 
statistical and elementary particle physics [201121) . The 
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Fig. 11. Finite-size scaling function U{tL^''^;p) of the ratio Fig. 13. p-dependence of the inverse of critical exponent 1/u. 
S2/S1 at p = 1.0. The inverse of the critical exponent l/v = The Monte Carlo data at p = 0.5, 0.6, 0.7, 0.1 
0.93 is taken for the scaling variable. are shown. 
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Fig. 12. Finite-size scaling functions U{tL^''^;p) of the ratio 
S2/S1 atp = 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. The corresponding 
1/1/ for the scaling variable at different p are presented in Fig |13l 
and Tabein 

Table 1. Critical reduced edge number re, ration of critical 
exponent jS/v and inverse of critical exponent l/f at different 
probability parameters. We obtain the critical reduced edge 



.(1) 



number rj, from S2/S1 and rj; from Insi 



(2) 



.(1) 



.(2) 



13/u 



1/u 



0.5 0.5000(2) 0.5000(4) 

0.6 0.5082(2) 0.5082(4) 

0.7 0.5153(3) 0.5153(2) 

0.8 0.5208(1) 0.5207(2) 

0.9 0.5244(1) 0.5244(1) 



0.108(4) 0.75 

0.102(3) 0.77(1) 

0.092(4) 0.79(1) 

0.081(5) 0.83(1) 

0.070(3) 0.88(1) 



1.0 0.5266(1) 0.52655(4) 0.064(3) 0.93(1) 
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Fig. 14. p-dependence of the finite-size scaling function 
U{tL^''^;p) of the ratio S2/S1 at the critical point i = 0. 



universality is characterized by the dimensionality d of 
the system and by the number n of the components of 
the order parameter j22| . Within a certain (d, n) univer- 
sality class, the critical exponents are independent of mi- 
croscopic details and are universal. In a finite-size system 
near its critical point, there is also universality. For ex- 
ample, the Binder cumulant ratio of magnetization at the 
critical point is universal. The ratio S2/S1 here is similar 
to the Binder cumulant ratio of magnetization. We could 
suppose that the ratio S2/S1 at the critical point is also 
universal and does not depend on microscopic details. 

In our previous investigations of the two-dimensional 
lattice networks under a generalized Achlioptas process, 
the dimensionality d of systems is fixed and the macro- 
scopic property of their order parameter is unchanged. 
But we have found in Figs l9ll3l and [U that the critical 
exponents /3, ly and the ratio S2/S1 at Tc depend on the 
probability parameter p. So the universality of percola- 
tion phase transition in these networks is characterized in 
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addition by the probability parameter p of GAP. A dif- 
ferent probability parameter p of GAP generates a differ- 
ent probability distribution of configuration. So the prob- 
ability parameter p is actually related to the macroscopic 
character of network. For a general classification of uni- 
versality class in complex networks, further investigations 
are needed. 



6 Conclusions 

We have investigated the percolation phase transitions 
in two-dimensional lattice network under a generalized 
Achlioptas process. In this GAP, we choose randomly two 
unoccupied edges in a two-dimensional lattice and the 
edge that minimizes the product of the two connecting 
cluster sizes is taken with a probability p. Our model 
becomes the two-dimensional bond percolation model at 
p = 0.5 and the two-dimensional lattice network under the 
minority product rule at p = 1. 

The size Si of the largest cluster in the lattice increases 
with the edge number Nr. When the reduced edge number 
r = Nr/N is larger than a certain value Tc, 5*1 becomes 
comparable with the lattice size N = L"^. At Tc, a gi- 
ant component emerges and there is a percolation phase 
transition. From the finite-size scaling analysis of Si and 
the ratio S2/S1, we can conclude that this percolation 
phase transition is continuous at probability parameter 
0.5 < p < 1. The critical exponent ratio /3/i^ can be deter- 
mined from the power- law behavior of ^i at re- To obtain 
the finite-size scaling function of the ratio S2/S1 from the 
Monte carlo simulation data of different L with the scal- 
ing variable tL^^'^ , the critical exponent of correlation 
length 1/ can be fixed. We find that the critical reduced 
edge number Vc increases with the probability parameter 
p, which is shown in FiglS] The critical exponent ratio P/v 
and the critical exponent i' decrease with the probability 
parameter p, as demonstrated in FiglHlandlTSl Under the 
GAP with 0.5 < p < 1, the formation of larger cluster is 
suppressed and this suppression increases with p. So the 
formation of a giant component should be delayed at a 
larger probability parameter p. This delay is accompanied 
then by the stronger emergence of the giant component, 
which is characterized by smaller /?. It is plausible that re 
increases and /3 decreases with p. The finite-size scaling 
functions of the ratio S2/S1 are given for different p in 

FiglH 

Within a certain universality class characterized by 
the dimensionality of the system and by the number of 
components of the order parameter, the universal quan- 
tities (critical exponents, amplitude ratios, and scaling 
functions) of different systems are identical. For the two- 
dimensional lattice networks under a GAP we discuss here, 
the critical exponents /?, v and the ratio S2/S1 at the criti- 
cal point depend on the the probability parameter p, which 
has been pointed out above. So the universality class of 
the percolation phase transition in this model should be 
characterized in addition by the probability parameter p. 
To understand the universality class of the critical phe- 
nomena in networks in general, further investigations are 



needed. For random networks, we introduce also a gen- 
eralized Achlioptas process and the phase transitions in 
these systems are investigated p5] . 
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